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We demonstrate by performing a global fit on Higgs signal strength data that large in- 
visible branching ratios (Br^) for a Standard Model (SM) Higgs particle are currently 
consistent with the experimental hints of a scalar resonance with mass ~ 124 GeV. For 
this mass, we find Brj nt) < 0.64 (95% CL) from a global fit to individual channel signal 
strengths supplied by ATLAS, CMS and the Tevatron collaborations. Novel tests that can 
be used to improve the prospects of experimentally discovering the existence of a Bii nv 
with future data are proposed. These tests are based on the combination of all visible chan- 
nel Higgs signal strengths, and allow us to examine the required reduction in experimental 
and theoretical errors in this data that would allow a more significantly bounded invisible 
branching ratio to be experimentally supported. We examine in some detail how our con- 
clusions and method are affected when a scalar resonance at this mass scale has couplings 
deviating from the SM ones. 

I. INTRODUCTION 

Two outstanding questions of importance that the LHC should shed light on are the origin of 
electroweak symmetry breaking (EWSB), and the relationship of the mechanism of EWSB to new 
states beyond the Standard Model (SM). 

There is strong indirect evidence for the EWSB sector being described by a theory that includes 
a particle that (at least) approximately has the properties of the SM Higgs boson. This evidence 
follows from many observables in flavour physics, from electroweak precision data (EWPD), LEP, 
the Tevatron and now the LHC. The SM Higgs is consistent with the results of these experimental 
probes in its pattern of breaking custodial symmetry (SU(2) C ) Q]-[3]| as well as in the manner 
by which it sources the experimentally established pattern of flavour violation. In light of these 
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results, there is substantial indirect evidence that a scalar field involved in EWSB will also be SM 
Higgs like in that the soft Higgs theorems of Refs. flU will be approximately respected, i.e. 
the scalar field will couple to the SM fields with a strength that is proportional to the mass of the 
corresponding SM particle. 

Directly, LEP, the Tevatron and LHC have jointly excluded large regions of possible Higgs 
masses in the SM. The upper bound (in the low mass region) for the SM Higgs is now restricted 
to m h < 130 GeV by ATLAS flU and m h < 129 GeV by CMS at 99 % CL, with a suggestive 
clustering of possible signal events around m h ~ 124 GeV. However, in spite of these results, 
the Higgs hypothesis is not yet established. In particular, there remains a significant freedom in 
the allowed couplings of a scalar effective field to the SM gauge bosons and fermions - so long as 
such a resonance has the approximate symmetries and properties discussed above [I8I4T01. 1 Such 
deviations in the properties of a scalar field from the SM Higgs can be interpreted as following 
from the Higgs boson emerging from a strongly interacting sector as a pseudo-Goldstone boson, 
or as the leading effect in the effective theory of more massive states that are integrated out. 

In light of this experimental situation, attempting to use current (and future) Higgs signal 
strength parameters to establish relationships between the EWSB sector and beyond the SM states 
is speculative. This is certainly the current status, as the suggestive clustering of signal excesses 
at nih ~ 124 GeV has not risen to the level of experimental evidence for a scalar resonance. 
Nevertheless, in this paper we will assume that future data will support the discovery of a scalar 
resonance at approximately this mass scale. Further, we will consider current signal strength mea- 
surements as indicative of the properties that such a scalar resonance has when performing global 
fits. In anticipation of such a discovery, it is of interest to consider how to efficiently extract 
evidence of yet other states coupled to such a scalar field. 

The gauge invariant mass operator of the scalar degrees of freedom, being of dimension two, is 
expected to couple generally to all degrees of freedom. It is difficult to forbid a coupling of this op- 
erator at the renormalizable (or non-renormalizable) level to new states. As such, the measurement 
of the decay width of a new scalar resonance to states that do not directly lead to significant ex- 
cesses in the Higgs discovery channels, defined in this paper as its invisible branching ratio Br in „, 
could be the first direct measurement of interactions with states beyond the SM. This exciting pos- 
sibility has lead to many studies on extracting the invisible width of the Higgs. In this paper we will 



For other model independent approaches to Higgs couplings determination see lfTTlfi"3ll . 
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explore a very straightforward route using current and future experimental results on Higgs prop- 
erties, expressed in terms of signal strength data, to probe for evidence of a Br inv . We first show 
(in Section [II]) that the current experimental hints of a new scalar resonance at rrih ~ 124 GeV do 
not put strong constraints on Bi inv . We then explore in detail how to extract evidence for (or ex- 
clude) a Bi inv using global combinations of best fit signal strength parameters, performing global 
X 2 fits, and demonstrating a global probability density function (PDF) approach that can be used 
to explore and optimize searches for Br inv in certain scenarios of beyond the SM (BSM) physics 



(Section [II A[ ). In Section [QBj we examine the related issue of the precision with which Bi inv is 
expected to be known in these scenarios when errors are small enough for a resonance discovery 
to be claimed. 

These promising results raise the question of the robustness of such a global approach. These 
techniques are most promising when BSM physics couples to the SM primarily through the 'Higgs 
portal' 03]-[T9]|, i.e. they are optimal in BSM scenarios where new states are not charged under the 
SM gauge group and couple to the SM (initially) through the SM gauge singlet scalar mass opera- 
tor. This is the case we explore in detail throughout Section [nj We briefly discuss and summarize 



the prospects for global fits to uncover Br int , in broader scenarios in Section III where the effective 
couplings of the Higgs to the SM fields deviate from their SM values due to the Higgs being a 
pseudo-Goldstone boson or due to the presence of higher dimensional operators. We present de- 
tailed numerics for these scenarios (based on current global Higgs fits to experimentally reported 



signal strengths) throughout the remainder of Section III 



We find that global fits to signal strength parameters will be a powerful approach to search 
for evidence of Bi inv in the scenarios we consider. But we note that challenges will exist in 
disentangling other new physics effects. This will likely require a combination of indirect global 
fit approaches to extracting Br inv , which is the focus of this paper, and more traditional direct 
searches for Br inv based on kinematic properties of Higgs signal channels. We compare and 



contrast these approaches in Section IV and then discuss our overall conclusion in Section |Vj 



4 



II. THE STANDARD MODEL HIGGS AND Br, 



tnv 



In this section, we will present the current global best fit 2 results for Br inv = T inv / (T inv + r SM ) 
for the SM Higgs, where T inv is the decay width to 'invisible' states, as defined above, and Tsm 
is the decay width of the SM Higgs. We perform a global fit to the available Higgs signal data, 
fitting to fifteen Higgs signal strength parameters /ij reported by ATLAS, CMS and the Tevatron 
collaborations which are defined as 



for a production of a Higgs that decays into the visible channel i. We use the best fit values of 
denoted by /tj, as reported by the experimental collaborations. The label j in the cross section, 
(jj^ h , is to denote that signal events in some final states are defined (by selection cuts) to only be 
summed over a subset of Higgs production processes j. We construct a global \ 2 measure on the 
fii by defining the matrix C as the covariance matrix of the observables, and A 9i = //j — jli as a 
vector of the difference between the signal strength variable m and the best fit value of the signal 
strengths 3 , 



Here i = 1 ■ • • N c h, where N c h denotes the number of channels. The matrix C is taken to be 
diagonal with the square of the 1 a theory and experimental errors added in quadrature for each 
observable, giving the error <Xj in the equation above. Correlation coefficients (currently not sup- 
plied by the experimental collaborations) are neglected. For the experimental errors we use ± 
symmetric 1 o errors on the reported /V For theory predictions of the Oj^ h and related errors, 
we use the numbers given on the webpage of the LHC Higgs Cross Section Working Group ll20l . 
The minimum (Xmm) ls determined, and the 68.2% (1 a), 95% (2 a) best fit confidence level (CL) 
regions are given by A% 2 < 1, 4, respectively, for x 2 = Xmin + ^X 2 - Here, the CL regions are 
defined by the cumulative distribution function (CDF) for a one-parameter fit. 

2 See Ref. [10] for a detailed discussion on our fitting procedure. The values of the fifteen input jli used are reported 
in the Appendix for completeness. Note that throughout this paper we will assume a 3% contamination due to gg 
events in the 77 j j signal strength, see the Appendix for further details. 

3 In a simple counting experiment the definition is jl = (n b s — nt)/nf M , in terms of the observed number of events 
(n i,s), the number of background events (n^) and the expected number of SM signal events (nf ). 



[Oj-+h X Br(/l i)} obse rved 

[<Jj^h x Br (ft, ->• i)]sM 



(1) 




i=i 



(2) 
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Although we are fitting for evidence of new 'invisible' states, we do not include effects due 
to new unknown interactions on the production of the SM Higgs in this section 4 . We include an 
invisible width by modifying the SM branching ratios universally for each decay into final states 
/via 



Thus, the effect of including an invisible width (of BSM origin) on the signal strengths is that 
the expected = 1 in the SM is modified to an expectation of //j = 1 — Br inv . We fit for the 
parameter Br inv assuming a SM Higgs with a total SM width T S m and a particular Higgs mass. The 
resulting \ 2 as a function of Bi inv is shown in Fig. [I] Interestingly, we find that the global x 2 is 
minimized for a non-zero value of Br inv . 5 An inspection of the data used for m h = 124 GeV (see 
the Appendix) reveals h — > 77 (ATLAS) and h — > WW (ATLAS and Tevatron) as the channels 
that most favor a nonzero value of Br inv , while h — >■ ZZ (ATLAS), h — >■ 77 (CMS) and h — > bb 
(Tevatron) are the channels which tend to drive Br inv — » 0. This result also demonstrates that, 
despite a suggestive hint in the data for a Higgs like scalar resonance, Br inv remains essentially 
unconstrained for the SM Higgs in the current data set. The allowed values are Bi inv < 0.37(0.64) 
at 68(95)% CL respectively for m h = 124 GeV, and Br inv < 0.39(0.65) for m h = 125 GeV. If 
the Br™ n > result (statistically marginal at this time) is confirmed by the future data set, and the 
Higgs is discovered, future global fits of this form could be the first evidence of the Higgs coupling 
to new states. 

It is instructive to look more closely at the x 2 fit to all channels we have performed rewriting 
Eq.[2]as 



In later sections, we examine the effect of contact interactions (that could be due to new states in a BSM sector) on 
our conclusions and method. In this section, we are essentially assuming the case of new states that are not charged 
under the SM group, coupling primarily to the SM through the SM gauge singlet operator H. 
5 Comparison of our results with those of Ref. ETI . which finds that a related global \ 2 is minimized for Br inv < 0, is 
not straightforward. For the results presented in this paper we use the fii (and errors) as reported by the experimental 
collaborations. The results in Ref. BTTl are based on /tj constructed from reported and expected CL limits that only 
approximate the experimentally reported /i,, apparently introducing a distortion in the data (and associated errors) 
that affect the final conclusions. Of course, due to the large experimental errors at this time, the 95% CL range is 
wide in our results and in Ref. ETI . 



Br{h -> /) 




(1 - Bi inv ) x Br SM (/i /). 



(3) 




(4) 
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FIG. 1: Global fit to the best-fit signal strength parameters in SM Higgs searches as supplied by the four 
experiments for individual channels (left) and for their combinations (right). These results are based on 
post Moriond 2012 data ( see the Appendix and Table I) when an invisible width is added to the SM as a free 
parameter. Left: The red (upper) solid curve isforrrih = 125 GeV; the blue (lower) solid curve isforrrih = 
124 GeV. The one sigma region defined with the CDF for a one parameter fit is given by the horizontal 
dashed lines in each case and the best fit points are given by (mh,^>Tinv) = (124,0.12), (125,0.15). Right: 
The red (lower) solid curve is for mh = 125 GeV; the blue ( upper) solid curve is for rrih = 124 GeV. Now 
the best fit points are given by (mh, Br^) = (124, 0.10), (125, 0.06). Comparing these results gives a 
sense of the effect of neglected correlations amongst the individual signal channels in such fits. 



where we have introduced the combined variables 



j N ch 
^2 = Yl ^2' 



N ch „ 



(5) 



Note that Eq. [4] is valid if all the ^ are equal, as is the case for the SM with an addition of Br ini) . 
This decomposition illuminates what our analysis of the fit to individual channels really does. 
The location of the minimum of the fit, and the Na intervals, is controlled by the first term in 
Eq. |4} which depends only on the combined parameters fi c and a c but not on the dispersion of 
the different /Vs around their average fi c . How good the fit is, is just given by the second piece 
in Eq.|4j which is simply Xmiw an d does depend on how separate are the individual channel /Vs 
from fi c . Interpreting, as we do in this paper, deviations of jl from its SM value of 1 as a Higgs 
invisible width, one immediately obtains that the x 2 is minimized (defining Br™ n ) when 



Br 



(6) 



This also offers the alternative approach of bypassing the individual channel analysis and using 
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directly the jl c values reported by the experiments in Table |lj We can use this data to do the \ 2 fit 
as the effect of Br im) on the signal strengths is a common multiplicative correction. 



Experiments 


fi c , m h = 124 


a c , m h = 124 


fi c , m h = 125 


a c , nih = 125 


CMS ma 


0.98 


0.30 


0.94 


0.32 


ATLAS El 


0.61 


0.38 


0.81 


0.38 


CDF&D0 EH 


1.31 


0.60 


1.28 


0.62 



TABLE I: Combined signal strengths and errors from ATLAS, CMS and the Tevatron collaborations. Here 
we quote ± symmetric la errors. 

These results lead to the combined values fi c = 0.9, a c = 0.22 (fi c = 0.94, a c = 0.23) for m h = 
124(125) GeV. This directly translates into the best fit results Br™ n = 0.09±0.22 (0.06±0.23) for 
m h = 124 (125) GeV (with the 95% CL limit Br inv < 0.54 (0.52)). Comparing the best fit value 
with the results of our previous analysis in terms of the individual channels: £l c are 0.88 (0.85) for 
m h = 124, (125) GeV. The results of the two fits, to the individual and to the combined /Ts are 
consistent within the quoted errors, indicating that neglected correlation effects in the individual 
signal channel fits do not dramatically change the fit results we will show. The slight preference 
for a nonzero invisible width is driven by ATLAS data at this time. 

A. Global PDF approach to discovering Br^ nv in the SM 

The x 2 approach we have discussed can be justified on the basis of a more detailed analysis 
that makes use of the combination of the PDF's for all sensitive Higgs search channels. To the 
extent that these PDF's are well described by Gaussian distributions, both approaches are basi- 
cally equivalent in terms of the discovery reach afforded. The combined PDF approach however 
has a more direct physical interpretation and makes clearer the expected experimental sensitivity 
required to discover or bound Br inv . Moreover, this treatment is more powerful as it could also 
capture possible deviations from simple Gaussian shapes. On the other hand, the \ 2 fit is useful 
to determine if a reduction in Higgs signal yields in all channels is really universal and leads to 
a value of Xmin indicating a good fit, and is very convenient for taking into account the effect 
of imposing EWPD constraints. In this section, we will discuss global searches from the global 
PDF perspective, discussing the relationship between the errors in fii and the discovery reach (or 
exclusion prospects) for Br im; with such an approach. 
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As in the previous x 2 analysis, the key point to the power of global searches is the fact that, 6 
in the presence of Br imj , the expected measured values of the strength of the signal with respect 
to pure SM expectations are modified as = 1) — > = 1 — Br im) ). Due to this, one can 
construct a global PDF (combining the PDF's of individual channels) sensitive to this shift in 
the signal strengths which aids in experimentally distinguishing Br in ^ ^ from the SM case 
Br ini , = . As in Ref. flU, we will assume that the PDF for each iii reported by the experimental 
collaborations can be approximated by Gaussian distributions 

pdf^^a^^e-^ 2 ^, (7) 

with one sigma error cr^, and best fit value fii for the signal channels i. This is the case (especially 
for fx near fii) as long as the number of events is large, > O(10) events, [|9l and systematic er- 
rors are subdominant. Using experimentally reported fii and <jj is the best and simplest way of 
approximating the likelihoods in the neighbourhood of fii. Experimental information on the 95% 
CL exclusion limits on iii give additional information on the PDF's up to higher values of ji and 
allow one to identify channels for which there are non-Gaussian tails. In such channels, using ex- 
clusion limits to extract the fii (as done in Ref. fl2TJ) will tend to overestimate the fii. This would 
subsequently bias an extracted value of Bi inv based on such constructed fii. 

A global combination of all the visible channels PDF's, where every PDF is approximated as 
above, is obtained as a product of the pdfi(Li) (where i = 1 • ■ • N c h) and it is also approximately 
Gaussian 

N ch 

pdf(Li, Ac, <Tc) oc J] pdfiQi, fii, a t ) = j\f c e -(«) 2 /M (8) 



where a c and fi c are defined in Eq.[5]and 



K V™ c 2 1 + Erf(/2 C /V2<7 c )' (9) 

Here we have normalized with the condition J °° pdf(n, fi c , cr c )dLi = 1 and Erf is the standard 
error function. In the limit where a, t m a, and correlations are neglected, one has the simple 



approximation a c ~ a/y/N ch . By combining all of the N ch visible channels, the distinguishability 
of Brj n „ > from alternative hypotheses (like pure background, or pure SM) is improved due to 
this ~ 1/ \/N ch suppression of a c compared to an individual visible signal channel's o"j. Assuming 



6 In the absence of an indirect impact of new states on production or decay of the Higgs through induced operators. 
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FIG. 2: Illustration of the probability density functions for the background-only (blue curve on the left) and 
SM ( red curve on the right) and corresponding p-values for a hypothetical fi c = 0.6. 

that with sufficient data the measured jl c converges to the theoretically expected fi c = 1 — Bi inv , 
by constructing a combined PDF one can determine the value of a c required to have the possibility 
to statistically pinpoint the presence of Bi inv ^ 0. 

To find evidence of a nonzero Br inv one has to be able to distinguish the Bi inv > hypothesis 
(dubbed SMinv) from the SM Higgs hypothesis with Bi inv = 0, and discern this case from the 
background-only hypothesis. We use the same approach used routinely in experimental analyses 
to estimate the significance of a signal excess in the data, which quantifies how unlikely such an 
excess would be if interpreted as an upward fluctuation of the background. One defines a p- value 
for the background-only hypothesis as 

POO 

Pback= pdf b {fi,a c )dfi , (10) 

where the background probability density function (or likelihood) is approximately a Gaussian 
centered at fx = (as n Q b s = n&) with some globally combined 1 -standard deviation spread (a c ) 
that results from the combination of N ch different channels each with an individual cr^: 7 

N ch 

pdf b ^,cr c ) « \[e-^l^ « e-^l^ . (11) 

i=l 

7 The overall normalization of pdfb(fx, er) can be fixed by J °° pdfb([i, er) = 1/2, which follows from 
J_ pdfb(fJ,, <j) — 1. Note that negative fi in pdfb have a physical interpretation in terms of downward fluctua- 
tions of the background. 
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For the i th channel, with an expected number of background events n^i and an expected number 
of signal events (in the SM) nff one has <jj = Jr\i/n^Y (neglecting systematic effects). Again, 
in the limit where all a { are comparable, and correlations are neglected, one has the simple scaling 
a c ~ o~i/y/N c h. A pi> ac fc-value as small as that corresponding to a 5a fluctuation will be required to 
claim Higgs discovery 8 . See Fig.[2]for an illustration of the definition of the background p- value. 

Besides having a small enough pback for a Higgs discovery, to claim evidence for Br inv > 
we should be able to discard also the pure SM hypothesis (with Br inv = 0), as a downward 
fluctuation in the signal yield could be misinterpreted as Br inv > 0. We proceed exactly as before 
and construct the global PDF for the SM hypothesis as a Gaussian centred at the SM value ji — 1: 

N ch 

pdf S M^, o c ) oc n e"^ 1 ) 2 /^ 2 ) = Msu e-^ 1 ) 2 /^ 2 ) . (12) 
i=i 

Here AfsM lS implicitly defined by the condition 



pdfsM^crc) = 1- (13) 
Also, neglecting systematic effects, cr^ = yjn b<i + n s s ^f /n s s ^ 1 . For nff 1 <ti n bji this is the same a; 



as in Eq. 1 1 so that we use the same notation for both. We then compute the p- value associated 



with the pure SM hypothesis as 



Psm = / pdf SM (^,a c )djJi . (14) 



o 



See Fig. |2]for an illustration. 

Claiming evidence for Br int , > requires having simultaneously a small Pback an d Psm- In 
order to quantify this, notice that pback < Pn requires fi c > N a c while psu < Pn leads to 
1 — Ac > N o c . Using fi c = 1 — Br inv , we can obtain, as a function of Br ira) , how small a c (the 
precision in the measurement of fi c ) is required to be for a Na evidence of nonzero Higgs invisible 
width. For reference, combining fifteen currently reported signal strengths from ATLAS, CMS, 
CDF and D0 , while neglecting correlations, one finds cr c ~ 0.3 for m h = 124 GeV. We estimate 
the current values of a c per experiment as cr c atlas — 0.5, ct C) cms — 0.4 and cr Cj Tcvatron — 0.6 
by combining the individual channels reported in the Appendix experiment by experiment (while 
neglecting correlations) using Eq. [5} This compares well to the combined fi c reported by the 



The p-value corresponding to an Na fluctuation is pjy = [1 — Erf (TV/ v2)]/2. One has p 5 = 2.87 x 10 
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experimental collaborations in Table [IJ Figure |3j left plot, shows the p-values for both hypothesis 
for a Br inv measurement with certain precision a c , chosen for illustration at around its current value 
~ 0.3 and future values a c = 0.15 and 0.05. As expected, claiming evidence for Bi inv > will 
be easier for Br inv ~ 0.5 and will be facilitated by a reduction in a c . With the current value, the 
plot also shows that the weak indication of Bi inv ~ 0.12 is perfectly compatible with a downward 
fluctuation of a SM-like Higgs, and even compatible with a background upward fluctuation at 
~ 3a. Figure [3](right), shows the required precision a c for la to 5a evidence of nonzero invisible 
width. As expected, the ability to find evidence of a nonzero Br inv degrades for small values 
of this parameter, when it is harder to disentangle SMinv from the SM and also for Br inv — > 1, 
when it is hard to discern a small signal over background. 9 As = ^/n^/n^f , with more 
luminosity the statistical component of a c will scale down with ~ 1/a/£ (so that the plotted l/of 
increases linearly with C). As an example, assuming this scaling of a c , and that the best fit value 
of Bi inv = 0.12 obtained in the global fit is the true value, this indicates that accumulated signal 
events should be increased (compared to the current data set) by a factor of « 25(100) to reach 
2(4) a evidence of this Br inv . 

At the end of the current LHC run it is expected that the accumulated luminosity will be enough 
to reach the level required for a 5a SM Higgs discovery per experiment. This means that both 
o"c,atlas an d o" c ,cms will be down to ~ 0.2 or lower. (This expectation is consistent with recent 
public statements by CMS and ATLAS, see Ref. Il22l .) Taking such values for these quantities, 
and combining with the current cr C) Tevatron> we arrive at a c ~ 0.15 (half the current value) as a 
reasonable number to expect by the end of the year. 



B. Bounding Br^ nv in the SM 

With the measured overall ji c , known with some error a c , we can also set 95% CL limits 
on Br inv . For this purpose one can use the overall PDF (from the combination of all Higgs 
search channels) for the signal strength parameter, which we again approximate by a Gaussian 

9 We have numerically cross checked the relationship between sensitivity to BRi„„ and a c shown in the p-value 
results with another simple test based directly on the lack of overlap of global PDF's. Introducing a PDF for the 
background only scenario, a SM PDF, and a test theory PDF, simply insisting that the N sigma allowed [i in the 
test theory PDF lies outside of the N sigma allowed regions of the other two PDF's, one finds a similar sensitivity 
to what is indicated for a 2 N a evidence for a common BR irl „ in the p-value test. 
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FIG. 3: Left: p-values for SM and background-only hypotheses (negative and positive-slope lines respec- 
tively) vs V>ii nv for several values of the la error a c on jl: its current value a c = 0.3 (solid); half of it, 
a c = 0.15, expected to be reached at the end of this year (dashed); and a future value a c = 0.05 (dotted). 
The red dashed horizontal lines show the p-values corresponding to significances from la to 5a 's. Right: 
a c (plotted as l/a%) required to pinpoint a non-zero Brj ni) with a significance from 1 to 5 a (curves from 
lower gray to upper red). The precise condition imposed is that, for a given p c = 1 — Brj n „, both Pback an d 
PSM (with Brj n „ = 0) are smaller than those corresponding to fluctuations from la to 5a. The horizontal 
dashed lines show again the values a c ~ 0.3, 0.15 and 0.05. 

centred at £i c with standard deviation o c . As we will interpret fi c < 1 as coming from a nonzero 
Br im , we restrict now fx to the interval (0, 1) and normalize the combined PDF accordingly, i.e. 
Jo pdf(fi)dfj, = 1. Then we determine a 95% CL interval (uli, V»L2) around fi c such that 

I pdf(n)dLL = 0.95 , (15) 

imposing the condition that the interval is centred at fi c if llli > and fi L2 < 1- Otherwise one 
fixes lili — or fj,^ — 1. From this interval we derive a 95% CL allowed band for Br im) as 

1 - A*£2 < Br im , < 1 - fMLt . (16) 

One can also place a comparable bound in the context of a x 2 fit that is given by 

Max [I- jl c - 2a c , 0] < BR im; < Min [l-fi c + 2a c , 1] . (17) 

The sensitivity of this % 2 -based bound is expected to be equivalent to the sensitivity to BR inv in 
the PDF test in the Gaussian limit. Figure [4] shows the sensitivity band, as a function of £i c for the 
PDF test, for several values of its error o c : the current one (a c = 0.3); the combined error expected 
when both ATLAS and CMS accumulate enough data for a ha Higgs discovery per experiment 
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1.0 

FIG. 4: 95% CL exclusion limits for Brj nt , as a function of the observed jl c for several values of its error a c : 
the current one (a c = 0.3); the combined error that is estimated to be reached by both ATLAS and CMS at 
the end of the year (o~ c = 0.15), with the 95% CL excluded region shaded; and with a future a c = 0.05. A 
red vertical solid line indicates the current value ofTSii nv obtained for = 124 GeV. The 95% CL limit 
for nih = 124 GeV obtained directly from the X 1 fit (< 0.64) is consistent with the PDF test results shown. 

over this year (<r c = 0.15), with the 95% CL excluded region shaded; and with a future error value 
down to a c = 0.05. As expected, if fi c is small this requires a large invisible width and a lower 
limit on Bi inv can be set while, if £i c is closer to 1, then only an upper limit on Br inv can be derived. 
For intermediate values of fi c a "measurement" of Br in „ would be possible. The plot shows that 
the error in the determination of the true value of Bi inv (along the diagonal) is approximately 2cr c . 
Note that for Br im) — > or 1, the corresponding values of fi c themselves require smaller o c than for 
moderate values of fi c to reach a discovery. Thus if a discovery is actually made in these cases, any 
corresponding Br^ will be simultaneously more accurately known than for the fi c ~ 0.5 case. 

Finally, in Fig. [5] we translate the best-fit value fi c obtained by combining the fi c results of 
ATLAS, CMS and the Tevatron into a best-fit value for the Br inv (using fi c — 1 — Bi inv ) as a 
function of m^. 10 The left plot shows the current situation, with a best-fit Br inv > for the 
interesting Higgs mass range m h ~ 124 GeV, which is in any case perfectly compatible with a 
zero value. The larger values of Br^ for other Higgs masses are also not statistically significant 
because they correspond to values of either psu or pback not particularly small (for reference, 
Psm < P2a only above the lower dashed line, while pbckg < P2a only below the upper dashed line). 

See also Ref. [23 ] for a recent analysis (on older data) with a similar reinterpretation of the data as in Fig. |5l (left). 
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FIG. 5: Left: Current status of the experimental situation concerning Brj nt „ extracted from combining the 
fb c values reported by ATLAS, CMS and the Tevatron, and interpreting deviations from fi c = 1 as coming 
from an invisible Higgs width. Above the lower dashed line Psm < P2a! below the upper dashed line, 
Pback < Via, so that no strong evidence for a nonzero value ofYiii nv is possible at this time. Right: Same 
as left, in a hypothetical future situation ( the solid curve is obtained from the left figure data series by 
shifting the data series to larger values ofYiii nv and reducing the error) assuming a factor 5 improvement 
in the precision with which the combined fi c could be measured compared to current data. The dashed lines 
correspond now to p-values equal to p^ a , so that finding 5a evidence for a nonzero Brj n „ would be possible 
in the white region between both lines. 

This can change with higher energy/luminosity and the right plot shows a hypothetical future 
situation with nonzero Higgs invisible width after collection of more data (such that the current 
a c ~ 0.3 used in the left plot is reduced by a factor 5). Besides a hypothetical curve with the best 
value for Br inv , the plot also shows the regions of parameter space for which p S M and pback are 
below 5(7, illustrating how such an analysis could claim indirect evidence for Br im , ^ 0. 



III. ROBUSTNESS OF GLOBAL FITS TO EXTRACT Br im; 

In the previous section we have examined the prospects for bounding or discovering Br in „ for 
the Higgs in BSM scenarios where new physics primarily couples to the dimension two scalar 
mass operator. In this section, we will examine how robust these conclusions are when a scalar 
resonance that has only approximately SM Higgs properties is involved in EWSB. First we will 



consider in Section III A the case of a minimal effective chiral EW Lagrangian with a non-linear 



realization of SU(2)l x U(l) y and a light scalar resonance. This scenario is most easily interpreted 
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in composite Higgs scenarios and introduces parameters (a, c) for the unknown coupling of a 
scalar resonance to the gauge and fermion fields of the SM, with the SM case corresponding 
to (o = l,c = 1). (See Refs. ||25"I - |27T0 . The obvious problem one faces in this case is how 
to determine if a universal reduction in signal yields is due to a nonzero Br inv or to a uniform 
reduction of the Higgs couplings involved in the search channels (i.e., a = c < 1). 

The robustness of global fits for Br inv in the presence of unknown higher dimensional operators 
is also important to determine. In Section [TTJ. when considering the effects of Br inv on the SM 
Higgs, it was assumed that all new light states are essentially not charged under SU(3) C x SU(2) L x 
U(1)y, so that only Higgs decay (but not Higgs production) was affected, due to invisible decays to 
those states. When this assumption is relaxed, one can consider BSM scenarios in which, besides 
the light SM singlets leading to Br inv , heavier new states [charged under SU(3) C x SU(2)l x U(1)y] 
leave their trace in the low-energy effective theory through higher dimensional operators. In that 
case, it is important to examine the interplay of the effects of such higher dimensional operators 
on Higgs production and the effect of the light singlet states on Br inv in global fits. We study this 



question in Section HID 



A. Global Fits to a Non-SM scalar resonance and extracting Br^ nv 

In this section, we consider the more general case, consistent with the current data set, of a 
minimal effective chiral EW Lagrangian with a non-linear realization of SU(2) x U(l)y and a 
light scalar resonance, denoted as h. Such a theory includes the Goldstone bosons associated with 
the breaking of the weakly gauged SU(2) L x U(l)y (which is a subgroup of SU(2) L x SU(2) R ) 
and the SM field content. This theory is the minimal effective theory of a light scalar degree of 
freedom that can have the experimentally supported pattern of MFV flavour breaking as in the SM, 
and respect custodial symmetry - SU(2) C , while the W and Z are massive. (See Ref. Il28l for an 
analysis of the LHC data relaxing the SU(2) C assumption.) The Goldstone bosons are denoted by 
7r a , where a = 1, 2, 3, and are grouped as 

E{x) = e laa7Ta/v , (18) 

with v = 246 GeV and a a the Pauli matrices. The T,(x) field transforms linearly under SU(2)l x 
SU(2)r as — > L Ftf where L, R indicate the transformation on the left and right under 
SU(2) L and SU(2) R , respectively, while SU(2) C is the diagonal subgroup of SU(2) L x SU(2) R , 
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under which the scalar h transforms as a singlet. The leading terms in the derivative expansion of 
such a theory are given by I1251 - I271 |29l 



C 



1 



1 + 2a 



h 



(19) 



V2 



1 + c 



h 



Vij Ur 

yfj d 3 R 



+ h.c. + 



Here we have neglected potential terms that are not relevant for the fits we will perform. Fitting 
the current data in such a theory has been recently explored in the literature fiHOH. Note that the 
SM Higgs is a special case of this theory, and corresponds to a linear completion (h becomes part 
of a linear multiplet) of this non-linear sigma model, with a = c = 1 . 



B. Imposing EWPD 

It is useful to consider fitting to EWPD simultaneously with the global data to obtain a more 
constrained parameter space in this effective theory. 11 In EWPD analyses, the corrections to the 
gauge boson propagators in this effective Lagrangian can be expressed in terms of shifts of the 
oblique parameters S and T ||30"] - |3~2]| given by 

AS » diz^l log (!*) , AT » JG-«1 log C'f) . (20) 
6 7r V A / 8 7r cos^ &w V A / 

These equations are approximate in that the numerical coefficient is determined from the loga- 
rithmic large rrih dependence of S, T given in Ref. [31]. Here we have introduced a Euclidean 
momentum cut-off scale A, which approximately represents the mass of new states that are re- 
quired to cut-off the growth in the longitudinal gauge boson scattering. In a full calculation, with 
all degrees of freedom, the cut-off scale will cancel. The degree to which this Euclidean cut-off 
properly captures the UV regularization of these integrals by new states not included in the effec- 
tive theory is model dependent. We assume that the UV completion of the effective Lagrangian 
is such that directly treating this cut-off scale as a proxy for a heavy mass scale integrated out is 
valid, i.e. that further arbitrary parameters rescaling the cut-off scale terms need not be introduced. 



The cut-off scale is chosen to be A = 4nv/\\/l — a 2 \ for 1. 



11 In fact, we will find in subsequent sections that marginalization over multiple parameters in the three dimensional 
fit space including EWPD is required to obtain a residual x 2 distribution that is not flat. 
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For EWPD we use the results of the Gfitter collaboration ll33l for = 120 GeV, 

S = 0.04 ±0.10, T = 0.05 ±0.11, U = 0.08 ± 0.11, (21) 
and the correlation coefficient matrix is given by 

I 1 0.89 -0.45 ^\ 

C = 



0.89 1 -0.69 
y-0.45 -0.69 1 j 



(22) 



We shift these results to having the input = 124 GeV using the one-loop contribution of 
the SM Higgs field to S and T. This numerical shift is ^ 10~ 2 . There is a strong preference for 
a ~ 1 in a global fit due to EWPD, i.e the SM mechanism of mass generation of the W ± and Z 
is strongly preferred in minimal scenarios where EWPD can be directly interpreted to dictate the 
value of a. When EWPD is imposed one has a bias in the fit space so that a > 1, but this should 
not be over-interpreted. This bias could in principle be a hint for the existence of other states 
in EWPD, but this possibility cannot be disentangled from cut-off scale effects without further 
experimental and theoretical input. We conservatively consider this bias to be simply a numerical 
artifact of our cut-off procedure. 

Interestingly, EWPD offers a handle to disentangling the degeneracy between a = c < 1 and 
the presence of Br inv : a ^ 1 has a direct impact on EWPD, while the new singlet states (into 
which the Higgs can decay invisibly) can have no impact on EWPD. In any case, the possibility 
of such degeneracy implies that further cross-checks of Bi inv > would be needed to confirm 
an eventual indirect evidence coming from the global tests discussed in the last Section. Directly 
confirming such indirect evidence for Bi inv is best accomplished in more traditional studies of 
experimental sensitivity to Br inv based on the kinematics of Higgs decay products. We discuss 



prospects for such a direct confirmation in Section IV 



C. Marginalizing/Fixing Parameters 

First, consider the case of fixing or marginalizing over one of the parameters (a, c, Br inv ) n to 
examine the robustness of our global fit results when Bi inv = is assumed, as in Ref. lfi~0l . Fits 



In the remainder of the paper we will always choose the value = 124 GeV as we have shown that the fit 
results are not strongly dependent (considering current errors) on the chosen mass (when varied in the range 124 — 
126 GeV). 
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Global Fit + EWPD, Br inv = 



Global Fit + EWPD, Br,„ v = 0.12 
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FIG. 6: Global fit to best-fit signal strengths in SM Higgs searches, for various values ofYiii nv . In the upper 
plots EWPD is also imposed, while in the lower plots EWPD is not included in the global fit. Here the 
green region is the 65% CL region defined through the CDF for a two parameter fit. The yellow region is 
the 90% CL region and the grey region encloses the 99% CL region. Also shown as solid black lines are the 
95% exclusion regions (outside this line is excluded at 95% CL) in the parameter space using the procedure 
described in Appendix B ofRef. / I70l/ and the data in the Appendix. 

with various Br ini , as an input value are shown in Fig. [6j We find that, when EWPD constraints are 
incorporated into the fit, the c < minimum of the x 2 is preferred for larger values of Bi inv . This 
is easy to understand, as Br(h — > 77) depends on the interference of fermion and gauge boson 
loops with an interference term oc —a c. As the invisible width gets larger and the expected number 
of events in 77 final states decreases, negative values of c allow (by constructive interference) the 
number of 77 events to be larger and more consistent with the data which show an excess in a 
number of 77 subchannels. It is interesting that this is another example where the breaking of the 
approximate c -h- — c symmetry in the parameter space has a physical consequence. If (relative) 
excesses in the signal strengths of h — > 77 became statistically significant with a larger data set, 
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FIG. 7: Allowed parameter space of (a, c) f/e/tj comparing the fit with EWPD (dot-dashed and dashed 
lines for 65% and 90% CL contours) and without EWPD. The plot colour convention is the same as in 
previous figures with the solid black line again denoting the 95% CL exclusion limit. Shown on the right is 
the allowed parameter space in the scaling variables (a — Br^m,, c/ y/\ — Y$ii nv ) with a solid line from 
(1, 1) to (0, 0) to mark the location of the SMinv point as a function of3ri nv . The dots on the line represent 
Brj TO from to 1 in steps of 0.1 from right to left. 

coincident with common suppressions in the other discovery channels, such a pattern of deviations 
can be explained with a large Bi inv and a negative c in the effective theory, with constraints from 
EWPD. Conversely, when EWPD are not used, the allowed parameter space is shifted to larger 
values of a and c as Bi inv increases to (partially) cancel the suppression of events due to Bi inv ^ 0. 
More precisely, the lower plots shown in Fig.[6]have a simple scaling property corresponding to a 
dilatation from the origin in (a, c) space, relating the spaces in plot i to plot j as 



{ai,Ci)yJl-BT\ nv = (a^cj)^/ 



(23) 



The constraints from EWPD on the fit space can be more easily understood by directly compar- 
ing the fit spaces as shown in FigJT] (left). Similarly, the dilatation scaling of the best fit space 
is illustrated in FigJT] (right) where we plot the fit space as a function of the scaling variables 

(a/ VI - Br inv , c/ VI - Br 

inv ) • 

Now consider treating each one of the parameters (a, c, Bi inv ) as a nuisance parameter in turn. 
Doing so we can also examine the effects of an unknown parameter on the remaining fit space. 
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mh =124 GeV, Br inv marginalized 



c 




mh =124 GeV (& EWPD), Br inv marginalized 
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FIG. 8: Allowed parameter space of (a, c) in the global fit when Br^ is marginalized over, subject to the 
physical constraint < Brj ni) < 1. To the left EWPD is not included in the joint fit, for the right EWPD 
is also included in a global \ 2 - The blue solid line delimits the parameter space where the prior condition 
Bim„ = is satisfied (within the blue line). 

For example, in marginalizing over Bi inv we define a reduced x 2 function 

X 2 (a, c) = x 2 (a, c, Br inv (a, c)), (24) 

where Bi inv (a,c) is given by the solution of dx 2 (a, c, Bi inv ) / dBi inv = 0. Then the allowed 
parameter space is defined through the CDF for a two parameter fit, and we obtain the results in 
Fig. [8} Marginalizing over the parameters c and a we find the results shown in Fig. [9] and Fig. 
[10] respectively which demonstrate the correlation between the allowed Bi inv , and the allowed 
parameter space for the remaining unknown parameters. This correlation is due to the dilatation 



relationship shown in Eq. 23 



Finally one can marginalize two of the free parameters simultaneously in order to obtain the 
residual x 2 distribution to examine if the slight statistical preference for Bi inv > persists. In 
this case, one must impose EWPD to avoid a flat distribution in the remaining free parameter. We 



find the results shown in Fig. 1 1 Of most interest is the result of marginalizing over free gauge 
and fermion couplings, while imposing EWPD. In this case, one finds that the global fit in this 
theory is now Br inv = 0, with the 95% CL limit Br inv < 0.57, for a scalar mass of 124 GeV. 
Comparing this to the SMinv result of Section II, we see that an unknown a, c can remove the 
slight preference for Br inv = in the current global fits when the Xmin w ^ m c > is the global 
minimum. Conversely when the Xmin w ^ m c < is chosen, the slight preference for Br inv in the 
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m h =\2A GeV, c marginalized m h =124 GeV, c marginalized (& EWPD) 




1.15 1.20 



a a 

FIG. 9: Allowed parameter space of (a, Brj W ) in the global fit when c is marginalized over, subject to the 
constraint \c\ < 3.5. The solid blue line illustrates where this prior is saturated (above the blue line). Again, 
on the right (left) panel EWPD is (not) included in the global X 2 fit- The plot colour convention is the same 
as in previous figures with the solid black line again denoting the 95% CL exclusion limit in the parameter 
space, to compare with the best fit regions. 



m h =124 GeV, (no EWPD) a marginalized m A = 124 GeV, a marginalized (& EWPD) 

1.0 ; . , ■ ■ . _ u _l L0 | I 




c 



FIG. 10: Allowed parameter space of (c, Br^) in the global fit when a is marginalized over, subject to 
the constraint < a < 3. The solid blue line illustrates again where the prior is saturated: above the 
horizontal line, a = 3 while, to the left of the vertical solid blue line, a = 0. Note that, in the right figure, 
the prior is never saturated as EWPD forces a ~ 1. Again in the right (left) panel EWPD is (not) included 
in the global X 2 fit- The plot colour convention is the same as in previous figures, and the 95% CL exclusion 
contours are black solid lines. 
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m h =124 GeV, (c,Br inv ) marginalized m h =124 GeV, (a,Br inv ) marginalized 
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FIG. 1 1 : Marginalizing over the two dimensional space to obtain the x 2 distribution for a single parameter. 
EWPD is included in the global x 2 measure. The blue dashed horizontal line in each plot is the 68% CL 
(la) limit, the red dot-dashed line in each figure is the 95% CL (2a) limit. The CL regions are defined by 
the cumulative distribution function for a one parameter fit. The bottom figure shows the two curves from 
the nearly degenerate X 2 m in w ^ c > (unfaded lines for (x^ jn )ij or c < (faded lines for (Xmin)^)- The 
difference in the minima is negligible [(xLinh ~ (xLn)i]/(xLn)i = °- 01 with (Xmin)l slightly prefered. 
When the c < minima is chosen, and the marginalization is performed, the preference for Br i nv in the 
current data increases due to the interference effects previously discussed. 

data set is not removed. Then the best fit is Br in ^ = 0.21, with the 95% CL limit Br int , < 0.75. 
This result makes clear that performing such a two dimensional marginalization over (a, c) is an 
important cross check to see if any evidence of Bx inv ^ is robustly preferred in a future data set. 

D. Higher Dimensional Operators, Brj nu , and the SM Higgs. 

It is also of interest to consider the impact of possible BSM states on Higgs production and 
decay, when evidence for Br int) > emerges from global fits. The exact impact of BSM states on 
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Higgs phenomenology is model dependent. In this section, we consider the case where new states 
that are SM singlets lead to Br inv , and other new states, that are charged under SU(3) C x SU(2) L x 
U(1)y, or at least a subgroup of the SM group, lead to higher dimensional operators. Our aim is 
to examine the degree to which conclusions about Bi inv can be extracted from global fits in the 
context of unknown Wilson coefficients of the resulting higher dimensional operators. 

Assuming that these BSM states do not source CP violation, the operators of interest for Higgs 
phenomenology (in global fits to fii) are given by 

_ C WB9l92 Hi T a HB ^ W a^ (25) 

Note that g x , g 2 , gs are the weak hypercharge, SU(2) gauge and SU(3) gauge couplings and we are 
using the notation of Ref. 113411 . The scale A corresponds to the mass scale of the lightest new state 
that is integrated out. We are primarily interested in the effects on <J gg ^ h and as these are 

loop level processes in the SM, sensitive to BSM effects. As we expect loop level contributions 
to these operators from the BSM states, we rescale the Wilson coefficients as Cj = c.,/(167r 2 ) for 
j = G,W, B, WB and fit to combinations of Cj. 

Using the results of Ref. [34], the effect of these operators are 



°gg^h ~ ° gg^h 



1 ~ 

V Cn 

1 - (1.39 + 0. 10T 



A 2 



r _ V SM 



V 2 C-v 

1 + 0.15—^ 
A 2 



(26) 



Here c 7 = cw + cb — cwb and we have used m t = 172.5 GeV, rrih = 124 GeV and a s (172.5) = 
0.1095. The imaginary part of the numerical coefficients above comes from including the b quark 
loop correction (we use = 4.7 GeV) in normalizing the BSM effect to the SM amplitudes. 
Normalizing in this manner is done to reduce the SM dependence in the BSM correction when 
this rescaling is used in our fits, and a numerical value is used for Cg^. h . As in Ref. [34J, we 
have retained the two loop QCD correction to the SM matching of the h G^ u G AtJ,u operator in the 
m t — > oo limit in these numerical coefficients. Due to this choice, this correction cancels out (in 
the m t — > oo limit) of the overall coefficient of the BSM effects when multiplied by the numerical 
value of CTgg^h- This is a ~ 10% correction on the quoted numerical coefficient. Initial state 
radiation and vertex corrections to G A V G tlv are expected to be common multiplicative factors for 
the operator h G A u G AtJlU in the m t — > oo limit, and as such are not incorporated in the numerical 
factors multiplying cq, c 7 above. We will consider the parameter space where the SM is modified 
by these corrections and Br inv in this section, fitting to (v 2 c 7 /A 2 , v 2 c G /A 2 , Br inv ). The exact 
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FIG. 12: Global fits to cq, c 7 assuming a SM Higgs for Brj n „ fixed to (0, 0.12, 0.5). EWPD is not simul- 
taneously imposed. The convention for the plot regions is the same as previous figures with the black line 
delineating the 95% CL exclusion contour. In examining the allowed parameter space recall that the factor 
1 /167r 2 has been scaled out of the BSM contribution, so that large allowed values of v 2 c 7 /A 2 , v 2 cq/ A 2 , 
although difficult to model build, are still perturbative corrections to the SM. 

relationship between these parameters, if any, is model dependent and unknown. As such, we fit 
to the data assuming no relationship between the three parameters. 13 

The operators in Chd also affect Br(/t — > 7 Z), where a different combination of the Wil- 
son coefficients cw,cb,cwb enters. This branching ratio is subdominant to the Bi(h -)• 77) 
branching ratio. (Numerically the values are Br(h — > 77) = 2.29 x 1CT 3 , and Br(/i — > 7Z) = 
1.46 x 1CT 3 . However, recall that when looking for the Higgs, the Z decay has to be multiplied 
by Bi(Z — > 11).) We neglect these effects when fitting for the allowed parameter space. We 
also do not include the effects of these operators on h — > W W, Z Z as the SM contribution is 
tree level for these processes. Further, we also neglect effects due to higher dimensional operators 
possibly modifying the differential distributions of the Higgs decay products, indirectly affecting 
the fii through modifying the effective signal efficiency for specific kinematic cuts. Such effects 
are expected to be negligible compared to the current uncertainties. However, we do not neglect 
the rescaling effect on that is identified with the rescaling on o- gg ^ h in Eq. 26 We include 



this rescaling consistently, which has a non-negligible impact on all branching ratios through the 
modification of Tsm- 

We show in Fig{l2]the allowed parameter space when Br ira) is fixed to prior values of 0, 0.12 



13 See also Refs. JSj |2T] [33), for example, for recent fits to BSM higher dimensional operator Wilson coefficients 
based on Higgs signal strength parameters. 
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FIG. 13: Marginalizing over cq, c 7 in the global fit (assuming a SMHiggs) and constructing the residual x 2 
distribution for Brj ni , (left figure). We show the results of two marginalization procedures in this figure. The 
solid (unfaded) line corresponds to marginalizing over cq, c 7 without any prior condition on these Wilson 
coefficients imposed. Conversely, the solid faded line corresponds to marginalizing over cg,c^ with the 
prior that the absolute values of each of these Wilson coefficients is < 1. For comparison, the right figure 
shows the result of marginalizing over Brj ni , when fitting for cq, c 7 in the current global fit. The plot colour 
scheme conventions are the same as in previous figures. 



or 0.5. In Fig 13 the residual \ distribution for Bx inv is shown when cg,c 7 are marginalized 



over (left) and we also show the allowed &g, c 7 parameter space when Br inv is marginalized over 
(right) subject to the prior constraint < Bi inv < 1. These results show the significant impact of 
the higher dimensional operators, in scenarios consistent with the assumptions of this section, on 
attempts to extract Bi inv from global fits to Higgs signal strength data. 

Most notably, we find that the slight preference in the global x 2 distribution for Br inv > is 
removed when marginalizing over such unknown BSM effects in the current data set. This offers 
further caution to over interpreting the slight preference in the global x 2 distribution for Br inv > 
at this time. Although in Fig|12]the required Wilson coefficient c 7 to still obtain a good fit when 
Br^v ^> is large, we find that even restricting cq, c 7 to clearly perturbative couplings (< 1), 
expected in many models, the preference for a Br inv > is removed. This result is shown in 
Fig{T3](left). 
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IV. PROSPECTS FOR DIRECT CONFIRMATION OF Br inv . 

As stated above, there is a degeneracy between the case a = c < 1 and the existence of a 
non-zero invisible decay. The former leaves the branching ratios unchanged due to a common 
suppression factor in the couplings, while all production channels are suppressed by the same 
common factor. On the other hand, in the simple case of leaving the SM couplings unchanged but 
allowing for an invisible width, the production channels are unchanged and the branching ratios 
are affected as in Eq. |3} leading to a common overall suppression of production times branching 
ratio compared to the SM. If a signal is seen with suppressed event rate with respect to the SM 
expectation the degeneracy between these two cases can only be removed by observing directly a 
non-vanishing invisible decay. 

It has been shown that associated production with gauge bosons, weak boson fusion and as- 
sociated production with top quarks allows one to discover a Higgs boson decaying invisibly, 
and to probe the invisible branching ratio. The typical signature is large missing transverse en- 
ergy/momentum. Assuming an invisible branching ratio of 1, a Higgs boson with mass up to about 
150 GeV can be discovered in Higgs radiation from a Z boson at J £ = 10 fb -1 and y/s = 14 TeV 
Il36l443ll . At high luminosity this reach can be extended to ~ 250 GeV in associated production 
with a top quark pair [|43] - l45l . Weak boson fusion allows for the discovery up to 480 GeV with 
10 fb 1 integrated luminosity 11431 l46]448l . Assuming SM production, invisible branching ratios 
as low as 25% can be probed in weak boson fusion for a 120 GeV Higgs boson at J C = 30 fb -1 
and y/s = 14 TeV at 95% CL [|46l448l . In associated production with a Z boson, branching ratios 
down to 45% can be probed ll38ll4Tl l42l while associated Higgs production with a top quark pair 
probes invisible branching ratios down to 56% ll45l . 

Recent papers have investigated the potential of a 7 TeV collider in direct searches for an 
invisible Higgs boson lfT2l |49] - [52l . The invisible branching ratio of a 125 GeV Higgs boson 
produced in weak boson fusion with SM strength can be constrained down to ~ 40% at y/s = 7 
TeV and / £ = 20 ftr 1 (501. Monojet searches from CMS based on 4.7 fb" 1 [|53l constrain 
£ = a / gsmBR{U — > inv) down to 1.3 at 95% CL translating to the constrained value of the 
invisible branching ratio in case of SM couplings [j5TI . 

The claimed 95% CL limits on £ expected for ps 124 GeV from the direct searches dis- 



cussed are summarized in Fig. 14 In this figure we also show for a direct comparison the 2a 



sensitivity expected in the global test statistics we advocate for the end of this year, with a c = 0.15 
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m h ~ 124 GeV 




FIG. 14: The sensitivity of different analyses of direct invisible Higgs decays can be assessed by comparing 
their potential in setting a 95% CL upper limit on £ = a j ' o~sAlBR(h —> inv^j. In the SM, with only Jii{ nv 
added as a free parameter, one identifies £ = Br i nv . The vertical lines show this reach for the indicated 
analyses (see labels at the bottom of the figure ). Also shown is the expected 2 a sensitivity band for the SM 
( where again £ = Bi{ nv ) using the global PDF test we advocate to compare with the sensitivity of these 
searches. The sensitivity band corresponds to the p-value test shown in Fig^(right). The light blue shaded 
horizontal band is for a c = 0.15, expected at the end of this year, and dark blue shaded horizontal band for 
a c = 0.05. Blue vertical lines: y/s = 14 TeV, J L = 30/& -1 , in weak boson fusion, Zh and tih production 
(from Ref. I\43V ). Green vertical line: y/s = 7 TeV, J C = 20 fb~ l , in weak boson fusion [50]. Red vertical 
line: y/s = 7 TeV, f L = 20 fb~ l , in monojet searches / 177V . 



and a future value when a c = 0.05. One sees in this figure that the global test statistics are likely 
to offer a significantly improved reach for the existence of Bi inv in the data set after the 2012 run. 
However, to claim a discovery of Bi inv will require a combination of these global searches and 



direct kinematic searches, as we have demonstrated throughout Section III 



V. CONCLUSIONS. 



In this paper, we have systematically examined the potential of global fits to extract information 
on Bi inv in the present and future signal strength data sets of a scalar resonance. We have focused 
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on the case of a scalar resonance with mass ~ 124 GeV, and have performed a global fit to the 
SM Higgs using the current signal strength data set, demonstrating that current global fits find the 
95% CL limit Bi inv < 0.64 for m h = 124 GeV. We have also illustrated how these results are 
statistically limited at this time and that any statistically significant conclusion on best fit values 
of Br inv will require more data. We have developed a new approach to globally combining signal 
strengths using global PDF's to optimize searches for new states that couple to the SM through 
the 'Higgs portal'. These promising results have lead us to examine the ability of global fits 
to resolve information on Br inv in the presence of unknown new physics effects simultaneously 
impacting the properties of the Higgs. Although disentangling these effects would require further 
experimental input, our results make clear the correlations expected between interpreting a global 
fit as providing evidence of Br inv and the sensitivity of such claims to other (unknown) new physics 
effects in the scenarios we have considered. Although the current signal strength data set we have 
considered in our numerical investigations only offers marginal evidence for a scalar resonance 
(and its properties) with m h ~ 124 — 126 GeV, the correlations with other new physics effects and 
the tests for evidence of Br inv we have explored are of continued interest as the data set evolves. 



Appendix A: Data Used 

The data we have used in the global fits of this paper are summarized in the table below. 
Due to an apparent inconsistency in the ATLAS best fit signal strength plot for h — > b b and the 
corresponding ATLAS CL S limit plot (that is under investigation by ATLAS) we do not use the b b 
best fit signal strength value in the combined fit at this time. For the pp — > 7 7 jj signal of CMS 
we assume a 3% contamination due to gg Higgs production events so that the relevant signal rate 
is given by 

(0Ma gg ^ h + a jjh ) x Br(h -)• 77). (Al) 

Here ajj h is given by VBF Higgs production. We do not use sub classes of WW events due to 
the lack of experimentally reported contaminations of these signal strengths due to other Higgs 
production processes. Simultaneously using a global best fit value fi for 77 events (for example) 
while also using a best fit fi for a subclass of events, such as 7 7 jj can result in a double counting 
of signal strengths that would incorrectly bias the fit. We avoid such double counting in our use of 
CMS and ATLAS data as the photon classes we use are exclusive, but note that double counting of 
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this form is present in Ref. iBTI . making it difficult to compare results. In particular, to avoid in- 
troducing such a bias is why we use the experimentally reported global ATLAS /t 77 , as a complete 
set of subchannel di-photon signal strengths is not available (in contrast to CMS). This is also the 
reason that we do not simultaneously use constructed signal strengths /t 77 and /i 77 ,p T >40Gev (from 
fermiophobic ll54ll searches). These signal strengths are not independent mutually exclusive event 
classes, being derived from the same signal event data. Our approach to this issue is different than 
the approach of Ref. lETI . 



Channel [Exp] 


Al24 


Al25 




rm _^ w+ W~ [CDF<VD(/)1 R e f |[24ll 


0.35+J-gf 


°- 03+ 003 

— U.Uu 




fCDF&DOl Ref (241 
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1 4+1-3 

1 . i ± ri ft 
— U.o 


PP ~ 


+ w W* -+ ^+ i/ r P [ATLAS] Ref. E3 


o.i±g : ? 
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1 2+ 2 - 1 
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-»• Z Z* ^ £+ £- £+ f [CMS] Ref. 1551 
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6 +0 - 9 
u -°-0.6 


PP 


-^WW*-^t+vfP[CMS] Ref. ED 


U.D_q y 
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■> 77 [CMS], Cat.4/BDT3, Refs. (52 [60] 
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■> 77 [CMS], Cat.3/BDT2, Refs. E3 |60l 


o.o+ 2 - 9 


9 9+1-5 
z - z -1.4 


PP ~ 


> 77 [CMS], Cat.2//BDT1, Refs-HHiQl 


2 1 +L5 


6 +L0 

u - u -0.9 
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+ 77 [CMS], Cat.l/BDTO, Refs. US ED 


J-^-i.o 


9 1+2.0 
z - i -1.6 




-> 77 jj [CMS] Refs. E3|6l 


q 7+2.5 
°- ' -1.8 


o R+2-2 
D -°-1.6 



TABLE II: Summary table of reported best fit signal strengths for various Higgs mass values. 

For rrih = 125 GeV we use the public results presented at Moriond 2012 that split the 77 
signal events into four (multivariate boosted decision tree -BDT) classes that are not identical to 
the classes used for m h = 124 GeV. See the relevant experimental papers for the detailed class 
definition in each case, but note that the event classes are exclusive (though correlated) and can be 
combined directly in our x 2 procedure. Once again correlation coefficients are neglected as they 
are not supplied, but the effect of pseudo-correlations have been examined in Ref. |[T0l and the fit 
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m h = 124 GeV 



m h = 125 GeV 



WW 
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WW -t llvv 
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FIG. 15: Pictorial presentation of the data in Table |77|/or = 124, 125 GeV. The central value of each 
signal strength and corresponding error band (la) is shown. Also shown is the background only hypothesis 
as a vertical dashed line, the SM hypothesis as a vertical solid green line at jl = 1 and a shaded vertical 
band that corresponds to fi c ± a c as defined in Eq. [5] The values of (l c are (0.88, 0.85) for = (124, 125) 
and are shown as a vertical solid black line. These parameters are related to the best fit value of J5i{ nv 
through Br™^ 1 = 1 — £i c . The results shown here are consistent with the global fit. 

was found to be stable against randomly chosen correlations. Also we have found in Section [TT] 
consistent results between two different approaches to the fit of signal- strength parameters: using 
the individual channels or the combined results. This indicates that neglected correlations do not 
bias the fit results outside the quoted errors. 

Note that here we use a value of 0.35±^f for the pp ->• W + W~ [CDF&D0] result for 
mh = 124 GeV, unlike in Ref. ifTOl . where we used the same value as for mh = 125 GeV. This 
introduces a small interpolation error, but allows better agreement with global combined signal 
strengths reported by the Tevatron collaboration. 
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